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INTRODUCTION 


This is the seventh quarterly report on Contract NAS5-27382 
entitled, "Spectroradiome trie Calibration of the Thematic Mapper and the 
Multispectral Scanner System.' In this report, we summarize the 
reduction of the data measured on July 8, 1984 at White Sands, New 
Mexico. The radiances incident at the entrance pupil of the Landsat 5 
sensors have been computed for bands 1-4. When these are compared to the 
digital counts of the TM image, we will have the ground based calibration 
of this sensor. The image has been received from Goddard SFC and is 
presently being analyzed. 
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JULY 8, 1984 CALIBRATION 


Field Measurements 

This was our first opportunity to acquire on-site measurements at 
White Sands in conjuction with Thematic Mapper imagery from Landsat 5. 
Imagery was not available for previous trips due to cloud cover. 
Assisting in the field measurements were Harumi Aoki, Ken Castle, Barbara 
Capron, Madgeleine Dinguirard, Ren Holm, Ray Jackson, Carol Kastner, Amy 
Phillips, Rich Savage, and Phil Slater. The instrumentation on hand 
included two Barnes radiometers, the cart and yoke, an old model of 
Reagan's radiometer, both of the Castle spectropolarimeters, four 
polycorders, a printer, the Compaq computer, and two 2 X 2 ft standard 
reflectance panels. 

Sunrise on the morning of July 8tb was at 6:10 a.m. New Mexico 
was on Mountain Daylight Time (MDT) this time of year, as are the times 
quoted here. The Reagan instrument was set up and began acquiring solar 
irradiance measurements at 7:15 a.m. (airmass 4). Rich Savage took 
temperature, humidity, and pressure readings there at Chuck site. He also 
arranged for a nearby radiosonde ascent. These data are presented in 
Appendix B. 

Two test sites bad been laid out on the previous trip. Each was 
a 4 X 4 pixel grid, aligned with the east/west scan lines expected of 
Landsatr A road with a 120° bend separated the two sites, and will 
facilitate the identification of the sites on the digital TM imagery. Each 
site was measured with a Barnes radiometer. Starting at the center of 


each 30 X 30 meter pixel, five reflectance measurements were Taken of the 
gypsum sands, within an area of about 5. X 0.5 meters. Reflectance panel 
readings were taken periodically during the course of these measurements. 
(Both the BaSO,, and Halon panels were recalibrated by Che Nianzeng 
immediately upon our return.) The data were averaged and recorded on 
polycorders. The site to the North of the road was scanned using the 
Purdue radiometer (one of the Barnes), mounted on the cart with the BaS0 H 
panel. Data were taken from 11 to 1 1 ; 20 a.m. The south site was scanned 
simultaneously with an USDA Barnes attached to a yoke. A Halon panel was 
used as the reference here. In addition to the two 4 X A pixel areas, two 
small areas were scanned near the van. This was done between 10:20 and 
11:40 a.m. These areas were selected for their contrast, representing 
extremes of light and dark for the local area. They will later be used 
in conjunction with aerial photography to map the reflectance of the 
entire area. Finally, diffuse to direct measurements were taken between 
8:30 and 10:20 a.m. These were made by comparing the radiance reflected 
from the ground to that radiance measured when the sun was blocked with 
a styrofoam parasol. These data can be used for comparison with the 
radiative transfer code, to verify our atmospheric models. 

A helicopter overflight was arranged for this trip. Jack Rees and 
co-pilot Keys piloted the helicopter and Jason Penny, PFC, was the 
photographer. The flight was almost exactly one full hour. It had the 
dual purpose of recording radiance at intermediate altitudes, again for 
comparison to the radiative transfer code output, and photographing the 
sites. Several slides have been scanned, using a microdensitometer, and 
will be used to characterize the ground reflectance. Five rolls of 
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Ektacbrome, ASA 100, were shot. All photos were taken at 1 / 500 tb of a 
second, with a 200 mm focal length lens. A series were taken at 6,000, 
2,000, and 500 ft., AGL (above ground level). The photos were bracketed 
from F/ 16 to F/22, in half F/stop Increments. Those taken at 2,000 ft. 
were particulary suited to our needs. Eight colored ground cloths had 
been laid out to define the two 4XA pixel areas that were measured with 
the Barnes, We noted that the blue and orange ground cloths were the 
most visible from the air, as well as being the easiest to see on the 
color slides. 

The field tests of the newly constructed spectropolarimeters will 
have to wait for our next trip to White Sands. The helicopter instrument 
suffered from a power failure, and data from the ground instrument, used 
for solar irradiance measurements, were lost during the TRS-80 to Compaq 
transfer. Back-up equipment was, however, available that enabled us to 
conduct an approximate calibration described in the following. 

Langley Plot Computations 

The Reagan radiometer was cycled through its narrowband filter 
set 95 times during the course of the morning. Each data set included, to 
the nearest second, a start and finish time and a voltage reading for each 
of the 9 spectral filters. These data were used as input to a program 
which computed solar zenith angle from ephemermis data. Using a 
refraction correction, airmasses were then calculated for each 
measurement. Finally, a weighted least squares analysis was used to 
compute the slope of the natural log voltages versus airmass, at each 
wavelength. The computed slopes are equivalent to the spectral optical 
depths of the atmosphere, t ex t» at the time of overpass. Temporal 


stability of the atmosphere is assumed* These results are given in Table 
1. Also shown are the optical depth components tftie> T Ray> and t oz * The 
computation of these components is discussed next. 

Using measured atmospheric pressure (883 mbar), tRay can be 
easily computed. After subtracting this from x ext , a curve is plotted, as 
in Figure la, which contains only Mie and ozone components. To determine 
rMie a curve is fit through all the text*"TRay data points that do not 
include absorption. Normally this is done by submitting the optical depth 
values to a routine which fits the data to an equation of the form log 
TMie* a e + a i lc, g X+ a 2 (log X)2. As many of the spectral filter data sets 
were rejected, an alternate approach was taken here. These problems 
arose because the older of the two Reagan radiometers was used. On this 
instrument the beater was not functioning. The data are less reliable 
without the temperature stabilization, due to fluctuations in detector 
responsi tivity.. This problem is thought to have affected mainly the 0.872 
and 1.03 pm channels. 

Instead of our normal procedures, therefore, a manual fitting of 
the data was performed. A curve of the form log T^ie* 3 „ +a jl°g X was 
assumed (a straight line on this log-log plot). This is an approximation, 
valid only if the aerosols can be correctly modeled as obeying a Junge 
radial size distribution, dn/dr=c r'"( v+ l). In such a case the data would 
fall exactly on a line whose slope, a t , yields the Junge parameter v via 
the relationship a^-v+2. By using only the 0.440 and 0.780 pm data 
points, a slope wss determined. At these wavelengths the ozone 

absorption coefficients are small. Actually, the 0.440 and 0.872 pm pair 
is prefered, as the absorption coefficients are approximately equal. The 
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Table 1. Langley Plot Results. 


Data from Reagan Radiometer 
8 July, 1984 

Chuck Site, White Sands, New Mexico 


Latitude 
Longti tude 
Right Ascension 
Declination 
Difference (Dec) 
Earth-Sun Distance 
Pressure 


32.935° 

106.407° 

7.226228 hours 
22.365361° 
-419.3 arc-sec 
1.016701 AU 
883. mbar 


WAV(ym) 

T ext 

T Mie 

T Ray 

T oz 

0.4000 

0.4426 

0.0981 

0.3172 

0.0000 

0.4400 

0.3060 

0.0922 

0.2138 

0.0006 

0.5217 

0.1921 

0.0824 

0.1063 

0.0127 

0.6120 

0.1543 

0.0743 

0.0555 

0.0246 

0.6708 

0.1091 

0.0699 

0.0382 

0.0098 

0.7120 

0.1063 

0.0673 

0,0300 

0.0046 

0.7797 

0.0842 

0.0634 

0.0208 

0.0027 

0.8717 

0.0948 

0.0589 

0.0133 

0.0006 

1.0303 

0.1103 

0.0528 

0.0068 

0.0000 


log TMie - a 0 +a,log X in ym 
a, - -1.269 a t - -0.654 


29123.7 (n 2 — 1 ) 2 

TRay * p 

2949810 , 25540 

“ ' 6432 - b+ -U6^ + 4PF 
X in micrometers 


T oz>X “ N0Z*a x 

NOZ * columnar ozone 

. T oz ’°*‘ la fin . 21 3.2 

® b * c l 2 ym 

a » spectral absorption 


ma tm-cm 
coefficient 


Table 2. Spectral components for TM midband wavelengths. 


WAV (ym) 

T ext 

T Mie 

T Ray 

T oz 

E„(mW/cm 

0.486 

0.2339 

0.0864 

0.1421 

0.0055 

175.955 

0.571 

0.1745 

0.0777 

0.0735 

0.0232 

180.580 

0.661 

0.1226 

0.0706 

0*0406 

0.0114 

153.916 

0.838 

0.0773 

0.0605 

0.0156 

0.0013 

104.708 
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data at 0.872 pm for this date, however, is unreliable, due to the 
temperature problem mentioned above. With the given constraints, the 
slope was found to be aj«-0.654, thus v*2.65, and a fl — 1.269. With these 
constants the TMie» TR a y> a nd ^oz components can be computed for any 
wavelength. Figure lb shows how each of these components contributes to 
the total optical depth. Table 2 gives the respective components for the 
TM midband wavelengths. These data will be used as input to the 
radiative transfer code. To be complete, a component of TH20 should be 
included for band 4. The contribution due to water vapor has not as yet 
been assessed. Also reported is the exo-a tmospherlc solar irradiance, as 
determined from interpolating the daca of Neckel and Labs (1981). The 
midband wavelengths of the TM sensor was computed by Palmer (1984) using 
the moments method and the pre-flight filter transmittance data. 

Panel Calibration 

The Herman radiative transfer code requires that the absolute 
reflectance of the gypsum sands be known. We have instead chosen to use 
the reflectance factor R(0 Z /O°). (A discussion of reflectance 
nomenclature and definitions is given in Appendix A.) Here 0 Z is the angle 
incident upon the gyspum, and 0° is the reflected angle, equal to the 
Thematic Mapper nadir-look angle. By using this quantity, the amount of 
light reflected in the direction of the TM is accurately characterized. A 
full BRDF characterization would be prefered. The gain in accuracy is not 
warranted, however, as the BRDF data would be difficult and time 
consuming to obtain. 

In the field the reflectance factor is measured with one of two 
Barnes radiometers. As these are uncalibrated they must be used in 


conjunction with a reference panel. The reflectance factor of the sands is 
determined via the relationship 

v sand * R ref 


^sand " 


^ref 


( 1 ) 


Were V sanc i and V re f are the output voltages of the' Barnes when looking 
over the sands and reference panel, respectively. These voltages are 
proportional to the radiance scattered upward and within the instrument's 
15® field of view (FOV). R re f the reflectance factor of the panel, as 
determined in the laboratory. On July 8^ one of the radiometers was 
assigned the Halon panel; the other radiometer was assigned the BaS0 H 
panel. While looking at the sands, each Barnes was periodically swung 
over the reference panel and voltage readings were recorded. Upon our 
return Che Nianzeng calibrated both panels. 

The calibration of the panels was conducted at the Optical 
Sciences Center, in a manner illustrated in Figure 2. A tungsten lamp was 
put at the focal point of an off-axis parabolic mirror. The emerging 
planar wavefront thereafter illuminated the reference panel at a known 
angle. The radiance reflected in a direction normal to the surface was 
then measured using a radiometer built by Che. Thereafter the reference 
panel was removed, and a primary standard surface put in its place. This 
primary standard was a Halon panel which had been calibrated by NBS on 
February 8, 1984. (They determined the reflectance factor RnbS^VO 0 ) 
for this surface.) 


In the first phase of the field panel calibration, the reflectance 
factor Rref(45°/0 ®) was determined. This was computed from the following 


A 0+k A 4* « < a* 



Voltmeter 


Figure 2. Laboratory set-up for the calibration of field reflectance 
panels. 
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Figure 3. Schematic of Chuck Site test area. 
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Table 3. Laboratory calibration of BaSO* panel. 
THE REFLECTANCE- FACTOR OF BaSC^ PANEL N0,S 


Irradiance 
Angle (deg.) 

430-470nm 

530-570nm 

630-670nm 

83Q-870nm 

10 

1,0420 

1.0355 

1.0205 

0.9858 

15 

1.0196 

1.0138 

0.9990 

0.9651 

20 

1.0011 

0.9945 

0.9807 

0.9479 

25 

0.9836 

0,9779 

0.9642 

0.9331 

30 

0.9668 

0.9610 

0.9479 

0.9185 

55 

0.9499 

0.9449 

0.9327 

0.9045 

40 

0.9326 

0.9279 

0.9169 

0.8905 

45 

0.9155 

0.9119 

0.9007 . 

0.8773 

50 

0.8971 

0.8937 

o.&asi 

0.8627 

55 

0.-8783 

0.8765 

0.8676 

0.8474 

60 

0.8597 

0,8582 

0.8498 

0.. 8326 

1 62 

0.8380 

0.3390 

0.8310 

0.8157 

70 

0.8157 

0.8X74 

0.8127 

0.7986 

75 . 

0.7926 

0.7929 

0.7892 

0.7789- 


I 

I 

I i 

| Data: July 12.1984 

I Sample: BaSO^ Panel No, 5 

1 Reference: Halon calibrated at (45/0) geometry by NBS (Feb. 

I * 8,1984) 

Location: Infra. Lab.,0SC 
| Viewing Zenith Angle: 0 deg. 

I Irradiance Angle: 10-75 deg. 


1 



Table 4. Laboratory calibration of Halon panel, 
THE REFLECTANCE FACTOR OF HALON PANEL CRAY’S) 


Irradiance 

430-470nm 

530-570na 

630 -67 Own 

830-87Qnm 

Angle (deg) 





10 

0.9965 

1.0007 

1.0020 

1.0042 

15 

0.9892 

0.9949 

0.9951 

0.9973 

20 

0.9829 

0.9872 

0.9888 

0.9901 

25 

0.9741 

0.9794 

0.9804 

0.9816 

30 

0.9648 

0.9702 

0.9716 

0.9728 

35 

• 

0.9556 

0.9597 

0.9618 

0.9636 

40 

0.9442 

0.9489 

0.9499 

0.9515 

45 

0.9319 

0.9368 

0.9377 

0.9391 

SO 

0.9166 

0.9210 

0,9231 

0.9253 

55 

0.8991 

0.9042 

0.9061 

0.9073 

60 

0,8784 

0,8832 

0.8860 

0.8877 

65 

0.8539- 

0.8608 

0.8633 

0.8638 

70 

<1.8221 

0.8297 

0.8323 

0.8343 

75 

0.7805 

0.7872 

0*7919 

0.7944 


Data: July 12, 1984- 
Sample: Halon panel f Ray's), 

Reference: Halon calibrated at (45/ 0) geometry by NBS (Feb. 8, 

1984) 

Location: Infra. Lab.,QSC 
Viewing Zenith Angle: 0 deg. 

Irradiance. angle: 10-75 deg. 


mm 

% 


V 


R re f(45 9 /0°) 


V r ef<45»/0 9 ) * Rnbs(45 9 /0 9 ) 


V NBS (45 9 /0 9 ) 

Next, the desired reflectance factor was computed from 


( 2 ) 


Rref(e/0°) 


V re f(e/0 # ) * cos45 4 * R re f(45 B /0 fl ) 


( 3 ) 


V r ef(A5 # /0») * cose 
The above steps were repeated for four spectral bandpass filters, 
each 40 nm wide. The measurement uncertainty was estimated to be less 
than IX. A small error in the reflectance factor was also introduced 
due to the non-uniformity of the panels. In the field the radiance 
reflected from the panels was averaged over a larger area. This was due 
to the 15* FOV, as compared to the 1° FOV of the laboratory radiometer. 
The results of the panel calibration are shown in Table 3 for the BaSO^ 
panel, and in Table 4 for Ray Jackson’s Halon panel. 


Reflectance Data 

Both the Barnes data and the panel calibration data were used to 
determine the reflectance factors of each of the test sites at White 
Sands. It is noted that the field measurements were taken at several 
times throughout the morning. Those data taken closest to the time of 
the Landsat overpass, 11:07:40 MDT, were used in the radiance 
computations. 

Several interpolations had to be made on the laboratory 
calibration data. The solar zenith angles were first computed for those 
times at which a field measurement of the sands was taken. The panel 
reflectance factors of Tables 3 and 4 were next interpolated to find the 
corresponding reflectance factors at these angles, for the four spectral 
filters available on Che’s radiometer. These wavelengths differ from 
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chose of Che Barnes, therefore one more interpolation was necessary. The 
intermediate computations were used to compute the reflectances for the 
eight Barnes' wavelengths. (The Barnes' wavelengths correspond to those 
of the TM, for bands 1-4.) 

The reflectance factors of the gypsum are determined from the 
above computed panel reflectances via Equation (1) above. A summary of 
this data reduction is given in Table 5. Only the data in channels 1-4 
are used here. Reflectance factors are given for pixels 1-16. A 
schematic of the site, which identifies these pixels is shown in Figure 3. 

Radiative Transfer Computations 

The mean reflectance values, given in Che bold type of Table 5, 
were used as input to the radiative transfer code. Also input were the 
atmospheric components listed in Table 2. The usual model assumptions 
were made for the aerosols. The maximum, minimum, and incremental radial 
sizes were 5.02, 0.02, and 0.04 pm, respectively. A refractive index of 
1 .54—0 -0 li was assumed. The code was run for both a solar zenith angle 
of 25° and 35®. The output, given in Table 6, is normalized for an exo- 
atmospberic solar irradiance of 1. After interpolating the data for a 
solar angle of 29.22° (that corresponding to the time of the overpass) the 
output is multiplied by the appropiate irradiance value. These final 
values are given in Table 7. They are the radiances that were incident on 
the TM sensor, for the morning of July 8. 

Summary 

Because of instrumentation problems and difficulties in accurately 
mapping ground reflectance due to the helicopter photography not being 










Table 5. Absolute Reflectance of White Sands test sites 


North Site 

TIME CH(1) 


ABSOLUTE REFLECTANCE 
CH(7) 


Ba S04 

10:52 

.957 

.952 

.939 

.915 

.858 

.839 

.715 

11:02 

.964 

.959 

.946 

.921 

.863 

.844 

.718 

ROAD 

10:54 

.510 

.577 

.619 

.650 

.625 

.527 

.259 

PIXELS 

10:57 

1-4,8 

.503 

.573 

.616 

.648 

.620 

.517 

.246 

10:57 

.510 

.581 

.627 

.660 

.632 

.527 

.236 

10:58 

.515 

.584 

.629 

.661 

.631 

.525 

.244 

10:59 

.530 

.601 

.645 

.680 

.646 

.535 

.231 

11:01 

.504 

.571 

.617 

.652 

.629 

.526 

,231 

BAS04 

11:02 

.964 

.959 

.946 

.921 

.863 

.844 

.718 

11:08 

.968 

.963 

.950 

.924 

.866 

.846 

.719 

PIXELS 

11:04 

7-5 

.475 

.538 

.579 

.610 

.589 

.492 

.233 

11:05 

.488 

.552 

.594 

.626 

.601 

.497 

.230 

11:06 

.504 

.571 

.613 

.643 

.611 

.502 

.224 

BAS04 

11:08 

.968 

.963 

.950 

.924 

.866 

.846 

.719 

11:15 

.973 

.968 

.954 

.929 

.869 

.849 

.719 

PIXELS 

11:08 

9-12 

.488 

.555 

.599 

.630 

.605 

.501 

.230 

11:10 

.504 

.574 

.620 

.653 

.623 

.511 

.229 

11:11 

.512 

.580 

.624 

.657 

.627 

.519 

.228 

11:12 

.526 

.597 

.643 

.676 

.637 

.520 

.213 

BAS 04 

11.15 

.973 

.968 

,954 

.929 

.869 

.849 

.719 

11.20 

.976 

.971 

.958 

.932 

.871 

.851 

.720 

PIXELS 

11.17 

13-16 

.501 

.566 

.608 

.637 

.611 

.512 

.252 

11.18 

.533 

.602 

.644 

.674 

.640 

.533 

.225 

11.19 

.511 

.581 

.624 

.656 

.631 

.526 

.228 

11.20 

.513 

.583 

.626 

.656 

.621 

.508 

.215 

CUMULATIVE: 
MEAN .507 

.576 

.619 

.651 

.622 

.516 


SDEV 

.015 

.017 

.018 

.019 

.015 

.013 






Table 5. Continued 


Site 

TIME 

Ha Ion 
10:58 

CRO) 

.965 

• • • 

.968 

of: pou 

ABSOLUTE REFLECTANCE 
CH(7) 

.970 .971 .973 .974 .980 

11: 

09 

.969 

.973 

.974 

.975 

.977 

.978 

.983 

ISC 

10 

: SCAN (PIXELS 1-] 
59 .489 .551 

L6) 

.599 

.642 

.625 

.500 

.233 

10 

59 

.494 

.557 

.601 

.638 

.618 

.497 

.241 

11 

00 

.491 

.548 

.590 

.627 

.604 

.481 

.229 

11 

00 

.468 

.529 

.573 

.612 

.591 

,468 

.229 

11 

01 

.487 

.546 

.588 

.626 

.600 

.480 

.246 

11 

01 

.493 

.550 

.592 

.632 

.606 

.482 

.228 

11 

02 

.509 

.572 

.615 

.653 

.630 

.508 

.251 

11 

02 

.514 

.576 

.618 

.658 

.628 

.500 

.231 

11 

03 

.486 

.546 

.590 

.629 

.614 

.491 

.236 

11 

04 

.503 

.567 

.611 

.650 

.634 

.519 

.256 

11 

04 

.516 

.578 

.622 

.660 

,637 

.512 

.243 

11 

05 

.498 

.559 

.603 

,641 

.621 

.500 

.241 

11 

05 

.510 

.572 

.616 

.654 

.624 

.492 

.228 

11 

06 

.510 

.569 

.610 

.647 

.616 

.489 

.230 

11 

06 

.506 

.567 

.610 

.650 

.632 

.514 

.250 

11 

07 

.503 

.561 

.603 

.643 

.622 

.500 

.250 

MEAN 

.499 

.559 

.603 

.641 

.619 

.496 

.239 

SDEV 

.013 

.014 

.013 

.013 

.013 

.014 

.010 

Ha Ion 
11:16 

.972 

.975 

.976 

.977 

.980 

.980 

.986 

11 

.27 

.976 

.979 

.980 

.981 

.984 

.984 

.989 

2N1 

11 

D SCAN (PIXELS 1- 
16 .492 .552 

16) 

.594 

.634 

.619 

.499 

.247 

11 

17 

.495 

.561 

.607 

.646 

.625 

.500 

.240 

11 

17 

.494 

.557 

.601 

.638 

.616 

.488 

.236 

11 

18 

.493 

.555 

.601 

.640 

.621 

.494 

.239 

11 

18 

.499 

.558 

.600 

.638 

.612 

.487 

.248 

11 

19 

.512 

.573 

.618 

.659 

.633 

.503 

.235 

11 

19 

.508 

.571 

.614 

.653 

.629 

.502 

.238 

11 

20 

.520 

.582 

.626 

.666 

.638 

.503 

.234 

11 

21 

.501 

.563 

.610 

.651 

.633 

.505 

.242 

11 

21 

.521 

.587 

.633 

.673 

.654 

.527 

.250 

11 

22 

.519 

.582 

.626 

.665 

.643 

.513 

.239 

11 

22 

.518 

.579 

.623 

.662 

.642 

.517 

.249 

11 

23 

.520 

.584 

.630 

.670 

.639 

.500 

.229 

11 

24 

.513 

.573 

.618 

.659 

.633 

.502 

.244 

11 

24 

.507 

.573 

.618 

.658 

.640 

.520 

.254 

11 

25 

.514 

.574 

.616 

.655 

.630 

.501 

.250 

MEAN 

.508 

.570 

.615 

.654 

.632 

.504 

.242 

SDEV 

.011 

.011 

.012 

.012 

.011 

.011 

.007 
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Table 5. Continued 


Site 

TIME CH(1) 


ABSOLUTE REFLECTANCE 
CH(7) 


BAS 04 

10:27 

.940 

.935 

.924 

.900 

.846 

,828 

.711 

10:32 

.943 

.939 

.927 

.903 

.849 

.831 

.711 

1ST DARK AREA 
10:28 .445 

.495 

.529 

.553 

.540 

.472 

.235 

10:28 

.459 

.508 

.541 

.562 

.543 

.473 

.228 

10:28 

.463 

.513 

.544 

.565 

.543 

.470 

.223 

10:29 

.458 

.509 

.54.3 

.564 

.543 

.470 

.224 

10:29 

.473 

.525 

.559 

.581 

.558 

.481 

.226 

WEAN 

.459 

.510 

.543 

.565 

.545 

*473 

.227 

SDEV 

.010 

.011 

.011 

.010 

.007 

.005 

.005 

1ST LIGHT AREA 
10:29 .516 

.580 

.623 

.656 

.607 

.480 

.178 

10:30 

.501 

.560 

.600 

.632 

.585 

.463 

.169 

10:30 

.504 

.564 

.602 

.635 

.589 

.469 

.174 

MEAN 

.507 

.568 

.608 

.641 

.594 

.471 

.174 

SDEV 

.008 

.010 

.013 

.013 

.012 

.009 

.004 

BAS04 

10:37 

.947 

.942 

.930 

.906 

.851 

.832 

.712 

10:40 

.949 

.944 

.932 

.908 

.852 

.834 

.712 

2ND DARK AREA 
10:37 .447 

.495 

.527 

.549 

.534 

.467 

.233 

10:38 

.455 

.505 

.538 

.559 

.543 

.473 

.233 

10:38 

.457 

.506 

.538 

.560 

.541 

.469 

.224 

10:38 

.457 

.507 

.540 

.562 

.541 

.467 

.224 

10:38 

.476 

.525 

.558 

.579 

.556 

.476 

.223 

MEAN 

.458 

.508 

.540 

.562 

.543 

.470 

.227 

SDEV 

.011 

.011 

.011 

.011 

.008 

.004 

.005 

2ND LIGHT AREA 
10:39 .516 

.578 

.618 

.651 

.604 

.477 

.177 

10:39 

.505 

.566 

.605 

.638 

.590 

.463 

.167 

10:40 

.504 

.561 

.597 

.628 

.583 

.464 

.172 

MEAN 

.509 

.568 

.607 

.639 

.592 

.468 

.172 

SDEV 

.007 

.009 

.011 

.011 

.010 

.008 

.005 

BAS04 

11:38 

.989 

.983 

.969 

.942 

.879 

.859 

.722 

11:42 

.991 

.986 

.972 

.945 

.881 

.860 

.723 
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W© '&*&' P 

OP POOH 


U 6 

QiUi k t 


3RD DARK AREA 
11:39 .478 

.534 

.568 

11:39 

.479 

.534 

.569 

11:40 

.487 

.543 

• 577 

11:40 

.485 

.542 

.577 

11:40 

.499 

.554 

.587 

MEAN 

.486 

.541 

.576 

SDEV 

.008 

.008 

.008 

3RD LIGHT AREA 
11.41 .540 

.608 

.651 

11.41 

.530 

.597 

.639 

11.42 

.530 

.598 

.639 

MEAN 

.533 

.601 

.643 

SDEV 

.006 

.006 

.007 


.589 

.569 

.495 

.245 

.591 

.567 

.488 

.233 

.598 

.573 

.493 

.236 

.598 

.571 

.490 

.233 

.608 

.580 

.495 

.232 

.597 

.572 

.492 

. *36 

.008 

.005 

.003 

.005 


.682 

.629 

.498 

.189 

.669 

.614 

.481 

.176 

.667 

.613 

.485 

.181 

.673 

.619 

.488 

,182 

.008 

.009 

.009 

.006 
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Table 6* Herman Code Output 


WAV 

TMie 

TRay 

tqz 

©z p 

E dir 

Edif 

Lpatb 

Lt 

0.486 

0.0864 

0.1421 

0.0055 

25 

0.507 

0.7001 

0.1519 

0.0246 

0.1334 






0.499 


0.1504 

0.0244 

0.1313 

0.486 

0.0864 

0.1421 

0.0055 

35 

0.507 

0.6156 

0.1417 

0.0223 

0.1190 






0.499 


0.1404 

0.0222 

0.1172 

0.571 

0.0777 

0.0735 

0.0232 

25 

0.576 

0.7477 

0.1113 

0.0134 

0.1457 






0.559 


0.1088 

0.0133 

0.1413 

0.571 

0.0777 

0.0735 

0.0232 

35 

0.576 

0.6621 

0.1028 

0.0121 

0.1299 






0.559 


0.1006 

0.0120 

0.1260 

0.661 

0.0706 

0.0406 

0.0114 

25 

0.619 

0.7916 

0.9078 

0.0080 

0.1618 






0.603 


0.8880 

0.0079 

0.1574 

0.661 

0.0706 

0.0406 

0.0114 

35 

0.619 

0.7053 

0.8337 

0.0072 

0.1447 






0.603 


0.8160 

C .007 1 

0.1407 

0.838 

0.0605 

0.0156 

0.0013 

25 

0.651 

0.8321 

0.6780 

0.0033 

0.1759 






0.641 


0.6683 

0.0032 

0.1730 

0.838 

0.0605 

0.0156 

0.0013 

35 

0.651 

0.7453 

0.6182 

0.0030 

0.1578 






0.641 


0.6095 

0.0029 

0.1552 


wbe re 



E dir 

The 

downward direct solar irradiance at the ground, cos6 z 

* exp(--r ex t sece z ). 

E dif 

The 

downward diffuse solar irradiance at the ground. 

Lpa cb 

The 

upward path radiance at the TM, Lx-CEtfir+Edif) * 
exp(-T exC sec5 # )p/ir 

Lt 

The 

total radiance at the TM a t a 5° nadir angle. 


Note all irradiance and radiance values are normalized for an exo- 
acmospberic solar irradiance of 1. 


Table 7. Computed radiance at Land, sat sensors 


L (mW/cm a sr ym) 


Band 

North Site 

South Site 

1 

23.18 

22.82 

2 

25.13 

24.37 

3 

23.82 

23.17 

4 

17.64 

17.35 


8 
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H 

il 
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vertical, the results of the July 8, 1984 measurement will have high 
uncertainties associated with them. However, the measurement attempt was 
worthwhile because of the experience gained in instrument operation and 
measurement. 

There are still a few remaining computations to be made with the 
July 8 C ^ data. Only a few discrete points of data were taken at the 
pixel centers of the two test sites. We are currently exploring the 
usage of the aerial photography (the slide imagery which has been scanned) 
to better characterize the reflectance of the area. The radiosonde and 
humidity data will also be looked at. Once the optical depth component 
T H20 can be determined, the Herman code will be rerun to account for 
atmospheric water vapor in band 4. The diffuse to direct data will be 
compared with the Herman code output. Finally, it remains to inspect the 
Landsat imagery, identify our test site, and compute incident radiance, as 
determined from the pre-flight data. A comparison can then be made of 
our in-flight calibration, to that made based upon pre-flight data and the 
internal calibrator data. The calibrations of bands 5 and 7 will not be 
investigated for this date. No atmospheric data are available at these 
wavelengths and -the optical depth components cannot be computed. The 
Castle spectropolarime ters will be equipped with a filter set which will 
allow us, in the future, to do calibrations in the short wave infrared. 
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APPENDIX A 


Reflectance 

The nomenclature, measuring geometry, and techniques associated 
with determining the reflectance of a given surface are quite varied. 
With this in mind, the technical basis of our measurements at White Sands 
is presented here. We are interested in the directional properties, as 
well as the magnitude, of the reflectance at the field site, for this 
reason we determine the reflectance factor. Reflectance factor and other 
quantities are defined below. A spectral dependence is assumed for each 
quantity. All of our reflectance measurements are made with finite 
spectral bandwidths. These are, in general, the 40 nm bandwidths of a 
laboratory radiometer, or those bandwidths associated with the Thematic 
Mapper and the Barnes modular multiband radiometer. 

Definitions and Nomenclature 

Reflectance Factor, (unitless). Ratio of the flux 

reflected by a sample surface to that which would be reflected into the 
same beam geometry by a lossless, lambertian surface which is identically 
irradiated. Thus, 

r 

Lt( 6' ,$')cos6' sine' d 0 ’ d$' 

, JlFOV , x 

R( 6 , 9 ' , q> ' ) * * (1) 

Lp( © 1 , <t> ' ) cos © 1 sine' de' d<$>' 

, IFOV 

is the reflectance factor measured with a detector having a given 


instantaneous field of view (JFOV). L$ is the radiance reflected off the 
sample target, and Lp is the radiance reflected off a perfect (lossless), 
diffuse surface. The incident beam originates from (e,$), and the 
reflected beam is viewed in the direction ( 0 1 , $') . 

Reflectance p (unitless). Ratio of the reflected flux to the 
incident flux. When referring to this parameter one needs to specify if 
the flux is integrated over the reflecting hemisphere, or if the reflected 
flux is measured within a given cone angle. The hemispherical reflectance 
can be related to the reflectance factor by 


2 ttJ 


it/ 2 


Lt( ®' j^Ocose' sine’ de' dV/E 


( 2 ) 


- R( 0 , 4>;2 tt) 

where E is the incident irradiance, generally from a well collimated 
beam. It is computed from 


f ( 


E - 


L c ( e,ij>)cose sine de d$ 


(3) 


1 J Ul 

Note the integration is over the solid angle doj*sine de d<f>. 

Bidirectional Reflectance Distribution Function (BRDF), f (sr”l). 
The ratio of the radiance reflected in the direction (e , ,<j> 1 ) to the total 
irradiance on the surface from the direction (6,$). 

f ( 6 , 4>; e ' , 4> * ) * Lt(6',$')/E (4) 

The quantity R(e,<j>;2ir) is equivalent to p. The 2ir denotes that 
the reflectance factor has been integrated over a hemisphere. R( 0 ; d) , or 
R(e/d), is an equivalent description, the "d” denoting that diffuse 
reflectance has been accounted for. Even when an integration is not 




l 


implied, the symbols <fr and are ofcen dropped for simplicity (as is true 
for any of the above parameters). 

Choice of the Reflectance Factor 

In calibrating the TM we are interested in knowing the radinnce 
reflected from the gypsum sands into a number of discrete angles. This 
allows both the directly and diffusely reflected solar radiation to be 
accurately characterized. Such a complete BRDF measurement is, however, 
both time consuming and difficult to measure. The equipment required is 
relatively complex, and there are difficulties associated with measuring 
the incident irradiance. Instead we have chosen to characterize the 
gypsum by the reflectance factor R(e z ;O fl ). This accurately describes the 
flux that is directly reflected towards the Landsat sensors. As the 
gypsum sands are not truely lambertian, some error is incured in not 
computing the full BRDF. Without the BRDF data, a lambertian surface is 
assumed. Thus an overestimate is made in the radiance not directly 
reflected towards the TM. Due to atmospheric multiple scattering in the 
atmosphere, some of this flux edventually reaches the sensors. This is 
the diffuse component of the radiance. The error made in predicting this 
term increases with increased multiple scattering”, and with departure 
from a lambertian behavior. Even so, the usage of the reflectance factor 
is justified. This is because the radiance received at the TM is dominated 
by the direct component, and multiple scattering is minimal for clear 
atmospheric conditions. 


Calibrating the Field Reference 

The reflectance factor Is measured with respect to a reference 

panel which Is calibrated In the laboratory to account for Its non-ideal 

characteristics. The calibration procedure was briefly described In the 

* 

body of the report. Here a development of the equations used in the two 
step calibration procedure Is given. 

To begin with, It Is assumed that a laboratory standard Is 
available. In our case a 50 mm diameter Halon disc was used which had 


been calibrated by NBS to determine R^bs^S 8 ;0°). The fictional parameter 
Vp(45 # ;0 8 ) Is thereby computed: 
v V NBS(45»;0») 

Vp(45 50 5 " RnBS<^ 6 ;0°) (5) 


This Is the voltage that would b :re been measured bad a perfect 
(lossless) lambertian surface been present. 

Using the above, the reflectance factor of l ie reference panel Is 
found for the same geometry: 


V re f(45° ;0 # ) 


( 6 ) 


VrefC^'iO 1 ) * RHjjsWiO*) 

' Vhbs<« •!(>•) 

In the next phase of calibration, the reflectance factor 
measurements are made at the angle of Interest, 6. For the Ideal 
lambertian surface the detector response at angle $ is easily predicted 
from the response at 45°. Such a surface reflects radiance uniformly 

into the upper hemisphere, thereby reflecting a factor of 1/ir of the 
incident irradiance. Thus, for this perfect (p-1) lambertian surface, 
illuminated with a beam of irradiance E(e), the following relationships 


flii tMMiiirT " f~1 i I ' mi 'in n ' pp . T - 






bold: 


Vp(45 # }0 # ) 

* R * E(45 J ) p/ir * R * E 0 cos45 # /ff 

(7) 


- R * E(6) p/t - R * E, cose/ ir 

(8) 


V p (45*;0 # ) * cose 



cos45 * 


The detector Is assumed to have a given response,, R, to the incoming 
radiance. This latter result is now substituted into the equation for 


RrefC^jO"), to yield the final, desired result: 


R ref(6;0°) 


Vref(9;0‘) 

Vp(e;0 9 ) 

Vref(0;O 6 ) * cos45 # 

Vp(45*;0*) * cose 

Vref(e;0 # ) * cos43 # * R re f<45 # ;0 # ) 
^ V re f(45 t ;0°) * cose 


(9) 


APPENDIX B 


Radiosonde and atmospheric data for July 8^, as provided by the 
Atmospheric Sciences Laboratory, White Sands Missle Range. 
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